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Abstract 



^ ; 

jy^ [ A nonrelativistic equation for the system of two interacting particles within 

the framework of a model with noncommuting operators of coordinates and 
momenta of different particles is proposed, and a self-consistent system of 

^ . equations for the wave function of every quantum state is deduced. Solutions 

K^ , for the lowest states of a hydrogenlike atom are found, and the comparison 

-^ ! with analogous solutions of the Klein-Gordon equation for the relativistic 

^^ '_ spinless problem is performed. In the case where the size of a two-particle 

system and the Compton wavelengths of particles forming it are of the same 

f^ . order, the essential differences with solutions of the Schrodinger nonrelativistic 

^^ . equation are revealed. 
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I. INTRODUCTION 

For a system of two interacting particles, the well-known Sclirodinger equation |l| can 
be derived through the formal substitution of the quantities E, ri, r2, pi, and p2 by the 
relevant operators on both sides of the relation E = i7(ri, r2, Pi, P2), where H is the 
classical Hamilton function of a system of two interacting particles. Here, ri and r2 are 
Cartesian coordinates of two particles with momenta pi and p2, respectively. 

The operators fi, f2, Pi, and p2 are such that they satisfy the following commutation 
relations: 

[xk, Pkx] = [Vk, Pky] = [zk, Pkz] = ift- {k = 1, 2) . (1) 

Here, h = h/27r, where h is the universal constant introduced by Planck. All other possible 
commutation relations equal zero, including 

[Xk, Plx] = [Vk, Ply] = [4, Plz] =0 (k^l) . (2) 

Equalities (^) are based on the assumption that measurements of coordinates and mo- 
menta of different particles do not disturb one another in principle even in the presence 
of some forces between particles 0. That is, one supposes that the change in the force 
action of a particle on another one, caused by a measurement of the coordinate of the first, 
propagates with finite velocity. 

Thus, to derive the Schrodinger nonrelativistic equation for a two-particle system, one 
uses, on the one hand, the Hamilton classical nonrelativistic function and, on the other 
hand, the implicit assumption about finiteness of the interaction propagation velocity. 

In the fully nonrelativistic quantum theory, we must consider the interaction propaga- 
tion velocity as infinitely large, which forces us to drop the requirement for the commutation 
relations (0) to hold. Having accepted this viewpoint, we will consider that, under a mea- 
surement of the coordinate of the first particle, there occurs the uncontrolled transfer of 
momentum not only to this particle but to the whole system since the particles are con- 
nected through the interaction potential, whose propagation velocity is infinitely large. 

Therefore, it is natural to require that the commutator of the coordinate operator of any 
particle and the of operator total momentum of the system Pc = pi + P2 be equal to ih: 

ih {k = 1,2). (3) 

Note that relations (^ hold true also for a Schrodinger nonrelativistic equation. Namely, 
they allow one to construct the operator of coordinates of the center of mass of the system, 
whose commutator with the operator of total momentum equals ih. On the contrary, the 
fulfillment of relations (|l]) is not obligatory for a system of interacting particles, and we 
intend to reject this requirement. 

The first attempt to construct the nonrelativistic equation for a system of interacting 
particles in the framework of the model with noncommuting operators of coordinates and 
momenta of different particles was undertaken in [^,|^ . However, the theory proposed there 
is phenomenological, i.e., it contains a parameter whose determination is ambiguous. More- 
over, the proposed self-consistent system of integro-differential equations was faced with 
difficulties of the probabilistic interpretation of a wave function. The present work is de- 
voted to the elimination of the drawbacks mentioned above. 
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II. FULLY NONRELATIVISTIC STATEMENT OF THE QUANTUM TWO-BODY 

PROBLEM 

The probabilistic interpretation of the square of the modulus of a wave function is pos- 
sible only under the assumption that measurements of coordinates or momenta of various 
particles do not principally disturb one another even if there exists some interaction between 
particles [Q. This means that the operators of coordinates and momenta of two particles 
commute with each other. 

As noted above, the operators of coordinates and momenta of different particles do not 
commute with one another in the general case. Let 

[xi, P2x] = ihfi , (4) 

where /i is a dimensionless Hermitian operator. Then it follows from Eq. (j^) that 

[xu Pix] = iHl - /i) • (5) 

By analogy, if 

[x2, Pix] = ifih , (6) 

then 

[X2, P2.] = tHl - h) ■ (7) 

The dimensionless Hermitian operators /i and /2 depend generally on the interaction force 
between particles F12 and on the masses of interacting particles mi and 7712- The operators 
/i and /2 cannot depend on a direction of the vector F12, since the commutation relations 
for the X, y, and z components should be identical analogously to (H)-(0), since there are no 
separated directions in the system, and the independent variables in the Cartesian coordinate 
system are fully equivalent. For this reason, the operators /i and /2 are only functions of 
the absolute value of a force, i.e., of |Fi2|: 

A = /l("^l,"^2, IF12I), /2 = /2(mi,^2,|Fi2|). (8) 

Let us make permutation of mi and 1712- Then 

[xi, P2x]=if^fi{^2,mi,\Fi2\) , [X2, Pix]=ifif2{m2,mi,\Fi2\) . (9) 

Compare (||) with (H), (H). Considering that the physical situation has not changed, we get 

/i(m2,mi,|Fi2|) = /2("^i,"^2,|Fi2|). (10) 

Thus, we have one unknown operator /i(mi,m2, IF12I). For 1712 — > 0, /i must tend to 
zero since, in the absence of the second particle, the whole momentum transferred under 
the measurement of the coordinate xi falls to the first one. If IF12I — >■ 0, then /i — ^ 0, i.e., 
without any interaction forces between particles, the operators of coordinates and momenta 
of different particles commute among themselves. The situation IF12I -^ 00 corresponds 
to the case where we have one particle of mass M formed by two strongly bound particles 



with masses mi and m2. Therefore, the momentum, received under a measurement of some 
coordinate, is distributed proportionally to masses of particles. This enables us to write 
down /i as /i = 7712/ M. Here, M = rrii + m2 is the system mass. 

Therefore, without loss of generality, we can present the operator /i as 

A = ^^(|Fi2|,mi,m2), (11) 

where £ is a new operator, which is assumed to be symmetric with respect to the masses 
of particles mi and 777,2. In what follows, we will omit its explicit dependence on masses to 
shorten formulas, namely, e(|Fi2| , m^i, 777.2) = ^(|Fi2|). For IF12I ^ 0, i ^ 0, and e — * 1 for 
IF12I -^ 00. 

For the noncommuting operators Xi and p2x, the uncertainty relation has the form 



A A ^ fim,2 
^ 2 M 
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(12) 



where (\1/ |e(|Fi2|)| ^P) ((^P | ^E'))^ = £ is a quantum-mechanical average of the operator i in 
the state ^^(ri, r2,t). 

Assume now that the two-particle system has only one quantum state \l/i(ri, r2,t). If 
we substitute the operator i in the formula for /i by its quantum-mechanical average 
(^1 |£(|Fi2|)| ^1) ((^1 I ^i))~^ = si, the uncertainty relation (|12D does not change. This 
allows us to construct a nonrelativistic equation for the two-particle system since the oper- 
ator /i is constant now. 

We present the commutation relations for all operators of coordinates and momenta in 
the two-body problem (for the y and z components, the commutators are the same) as 

[xi, pix] =ih(l- ^^1) ' (13) 

[x2, P2x] =ih(l- —eij , (14) 

[xi,P2x] = ih^eu (15) 

[X2, Plo;] = ^^^^1 ' (16) 

[xi,a;2]=0, (17) 

[pix,P2x]=0. (18) 

Now we can construct one of the possible representations for the operators of coordinates 
and momenta of the two-particle system as 

f 1 = ri , (19) 



Pi 



r2 = r2 



-^h[l- ^ei) Vi - ^h'^e,V^ , 



(20) 
(21) 



P2 = -«^^£i Vi -ih{l 
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(22) 



Here, as independent variables, we use the coordinates of the particles ri and r2 , since the 
corresponding operators commute with one another. 



By changing operators (p^-(p2D in the Hamilton function, we get the self-consistent 
system of integro-differential equations for the nonrelativistic two-body problem: 



d 
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where 



£1 
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if = f + I>(|r2-ri|), 



(24) 
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(26) 



We recall that we have received the self-consistent system of equations (p^)-Pq) for the 
two-particle system by assuming the existence of only one quantum state \I/i(ri, r2,t). 
Under the standard change of variables 



r = r2 - ri. 



(27) 



R 



miTi + in2V2 
mi + m2 



(28) 



where r and R are the relative distance between particles and the coordinates of the center 
of mass of the system, the system of equations (p3|)-(E6|) becomes simpler: 



z^-^i(r,R,t) = 



Ar ^ , ^' A, + y(rl 



2M 



2/i 



^i(r,R,t), 



(29) 



£1 



{^i\e{\Fi2{T)\)\^i) 
(^il^i) 



(30) 



Here, /i is the reduced mass, /i = 71111712/ {mi + 7712). The operator of total momentum reads 

Pc = Pi + P2 = -iWR. (31) 

Let the Hamiltonian H do not exphcitly depend on time. Then we obtain, by the 
substitution 



^1 = ^iexp ( -^^j 



(32) 



where ■j/'i depends on coordinates in the configuration space but not on time, the self- 
consistent system of integro-differential equations for a stationary state of the two-particle 
system as 



Ar ^ ^., ' Ar + V{r) 



2M 



2/i 



^i(r,R) =E^i(r,R), 



(33) 



£1 



(V;i|e(|Fi2(r)|)|^i) 
(^1 I A) 



(34) 



By the substitution ipii^r, R) =$(R) v^i(r), we can separate the motion of the center of 
mass of the system as a whole. As a result, we arrive at the following self-consistent system 
of equations: 



n\l-ei] 
2/i 



-Ar + V (r) 



V9i(r) =EiV9i(r), 



(35) 



£1 



(v9i|e(|Fi2(r)|)|v9i) 



(36) 



As in the Schrodinger nonrelativistic theory, a wave function ipi (r) should be continuous 
together with its partial derivatives of the first order in the whole space and, in addition, be 
a bounded single-valued function of its arguments. 

As in the Schrodinger theory, for particles interacting by means of a spherically symmetric 
potential, which depends only on the distance between particles, the wave function v'i(r) 
can be represented in the following form: 



1 /r 

<y5i(r) =-Xii{r)Yi^ 



r 



(37) 



where Yim{r/r) are orthonormalized spherical functions. Then the function XiK^) satisfies 
the following system of equations: 






2/i 



V{r) 



Xii{r) = EiiXiiir) 



(38) 
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(Xnkl|Fi2(r)|)|xn) 
iXii I Xii) 



(39) 



We emphasize that, for a given value of /, the corresponding value of eu can exist. 

Now we consider the laws of conservation which exist for the proposed fully nonrelativistic 
scheme of the system of two interacting particles. Since no external fields act on the system, 
its Hamiltonian must be invariant relative to a parallel translation of the coordinate system 
in space and time as well as relative to a rotation of the coordinate axes. In addition, the 
motion equations do not change under a uniform and rectilinear motion of the system (the 
Galilei invariance). 

The Hamiltonian of the isolated system p5| ) does not depend on time explicitly, and 
therefore the energy of the system is an integral of motion. 

The operator of total momentum of the system P^ = P1+P2 = —iTiS/i—iKV2 is connected 
with the operator of infinitesimal translation, which transforms the function \l'(ri,r2) into 
the function \E'(ri + (5r, V2 + 5r), by the relation 

2 i 

l + (5r-^Vfc = l + -<5r-Pe (40) 

fc=i 



and commutes with Hamiltonian (25): 



H.V 



0, (41) 



where 5v is the vector of infinitesimal parallel displacement of all radius- vectors by the same 
quantity: r^— i> r^ + 5y. Thus, three components of the total momentum are integrals of 
motion, and the total momentum of the two-particle system is conserved. 

By virtue of the isotropy of space, the Hamiltonian of a closed system must not vary 
under a rotation of the whole system by an arbitrary angle around an arbitrary axis. It 
is sufficient to require the fulfilment of this condition for any infinitesimal rotation, whose 
vector 50 is equal to the rotation angle 5VL in modulus and is directed along the rotation 
axis. The operator of infinitesimal rotation, which transforms the function ^(ri,r2) into 
\l/(ri + [(5f2 X ri], rg + [50 x r2]), is connected with the operator of total angular momentum 
of the system through the relation 

^ i 

l + 5n-Y. [vk X Vfc] = 1 + -5rj-L (42) 

fc=l '^ 

and commutes with Hamiltonian (^) of the system. Thus, the total angular momentum 
L = — ^^Z]fc=i[rA; X Vfc] of the two-particle system is conserved. 

It is worth noting the following. If we write the operators li = [fi x pi] and l2=[f2 x P2] 
in a formal way for each particle, it is easy to show that they and their sum are not angular 
momenta because they do not satisfy the standard commutation relations peculiar to angular 
momentum: 

[L^,Ly\=i%L^, [Ly,L^]=i%L^, [L^,L^] = ihLy . (43) 

However, from the operators f^ and p^, one can construct the operator possessing the 
mentioned properties of total angular momentum: 

2 
L = Zl Ckn[rk X p„] = -ihJ2[^k X Vfc]. (44) 

k,n fe=l 



For the two-particle system, the coefficients Ckn are as follows: 



Cu=[l-'^e,]{l-e,)-^ 



Cu = -"^6,(1 - e,r\ 



C,, = -^e,il-e,)-\ 



(45) 
(46) 
(47) 



G 



22 






ei 



(48) 



It is necessary to emphasize that the parameter of noncommutativity of the operators 
of coordinates and momenta of different particles ei is the quantum-mechanical average of 
an operator which depends on the modulus of the interaction force between two particles 
(i.e., on the distance between two particles). Therefore, this operator is independent of 
motion of the center of mass of the system since the system of equations ([23|)-(P^ admits 
the separation of motion of the center of mass of the two-particle system as a whole. 

Upon the motion of two reference systems relative to one another with constant velocity 
V, the operators r^ and p^ are transformed into f^ — vt and p^ — TTi^v. Such a Galilei 
transformation of the system of particles is described by the operator [Q 

(49) 



G(v,t) = exp zv- MR - PM /Ti 



where M, R, and Pc are the mass, operator of coordinate, and operator of momentum of 
the center of mass of the two-particle system. It is easy to show that the condition for the 
Galilei invariance of the equation of motion 



G-^(v,t) 



ih 



d_ 
di 



H 



G(v,t) 



ih 



d_ 
di 



H 



(50) 



is valid for Hamiltonian (PSI). 

The quantity eim2/M is the average share of the momentum transferred to the second 
particle under the measurement of the coordinate of the ffist. The value of this momentum 
can be estimated as |Fi2| At, where At is the duration of measurement of the coordinate of 
the ffist particle. Therefore, the operator e can be presented in the form 



^l|Fl2|) 



Fi.lAt 



12 



Fi2|At + Po 



12 



+ Fo- 



(51) 



where Pq or Fq is some constant of the relevant dimensional. For the two-particle system, we 
can construct the quantity Fq from two constants /i and h appearing in the problem under 
study as 



{^ic 



2\2 



he 



(52) 



i.e., this is the ratio of the energy of rest of the particle with reduced mass /i to its Compton 
wavelength. Therefore, the parameter of noncommutativity of the operators of coordinates 
and momenta of different particles can be presented in the form: 



£ii 



( 



V 



xli{r)- 



|Fi2(r) 



\ 



Fi2(r)| + 



{^ic 



,2\2 



■dr 



xli{r)dr 



(53) 



he 



J 



The proposed method of determination of the parameter of noncommutativity of the 
operators of coordinates and momenta of different particles contains no parameters and 
therefore substantially differs from the phenomenological method, proposed in [Q. In the 
latter method, the representation of the operator e = 1 — exp {—VLqF12 (|r2 — ri|)) includes 
the parameter VLq, whose determination is ambiguous. 

Now we consider the situation where a two-particle system has A^ quantum states, each 
of them is characterized by the wave function ^i, . . . , \I/Ar. Let the energy of every stationary 
state be denoted by £^i, . . . , Ejq, respectively. If we take the parameter of noncommutativity 
e to be the same for each state and equal to, for example, the maximum value of 



max 



(^l|£(|Fi2|)|^l) (^;v|^(|Fi2|)|^7v) 



(^1 I ^l) 



{^N I ^n) 



(54) 



then operators (pT|)-(p^) will be linear. In this case, the principle of superposition is valid 
for states. The probabilistic interpretation of wave functions is also conserved. 

However, such a choice of the parameter of noncommutativity e is not physically correct 
because there exist different average values of the modulus of the interaction force between 
particles in different quantum states, and the parameter of noncommutativity should be 
generally different in each quantum state. In this case, we have the own self-consistent 
system of equations of the type of (p9D-(pOD for each quantum state, i.e., each quantum state 
possesses the own interaction Hamiltonian. It is possible to reconcile this situation with the 
principle of superposition and probabilistic interpretation of wave functions analogously to 
the introduction of spin into the Schrodinger equation by Pauli . 

We recall that the motion of an electron in the constant homogeneous magnetic field Hz 
directed along the z axis is described by two equations with regard for the electron spin: 



H,+ 



eh 
2me 



H. 



^1 = ^1^1, 



(55) 



H. 



eh 
2me 



-H. 



^, 



Eo^ 



2, 



(56) 



where the wave functions \l/i and \l/2 describe the states with the projections of the spin 
onto the z axis s^ = +| and s^ = — |, respectively. Here, Hq stands for the Hamiltonian of 
the Schrodinger equation for a charged particle in the external electromagnetic field. The 
full wave function is two-component by the proposition by Pauli and is written in the form 
of a matrix column: 



* 



^1 

^2 



(57) 



As a function ^^, we choose the so-called Hermite-conjugate wave function ^^ = ( ^^, \E'2 
whose elements are not only complex conjugate but also transposed. Then the density of 
probability is given by 



^t^ = fx^*^ ^* 



^1 

^2 



\I/*\[/^ + \[^2^2- 



(58) 



Consider a system with A^ quantum states (A^ can be infinite, and some share of functions 
corresponds to a discrete spectrum whereas the rest to a continuous one). Similarly to (^7|), 
we introduce wave functions, each of them is represented by a matrix with one column and 
A^ rows, where \E'i, . . . , "^^ normed in a proper way describe possible states of the quantum 
system with the energies Ei, . . . , E^, respectively: 



*i 



/*i(r,R,t)\ 



, *2 = 


( ° ^ 

^2(r,R,t) 



, . . . , *7V = 


^ ° ] 


V ) 




I J 






^^jv(r,R,t)/ 



(59) 



In correspondence with such a construction of state vectors, the quantity ^k is related to 
that situation where the system has the energy Ej. with probability one. Then any quantum 
state of the two-particle system can be presented as 



* 



/ ai^i(r,R,t) \ 
a2^2(r,R,t) 



U7V^7v(r,R,t)/ 



(60) 



where a^ are arbitrary numbers, complex ones in the general case. The requirement that 
the state vector (|57D be unit, i.e., the scalar product of ^ and the corresponding vector ^^ 

o 

conjugated by Hermite be equal to 1, allows one to interpret \ak\ as the probabihty for the 
system to be in the state "^k'- 



I 



*^*(ir 



N 

fc=i 



1. 



(61) 



Now the wave equation for the fully nonrelativistic two-body problem has the following form: 



ih 



d_ 
di 



/vl>i(r,R,t) \ 
^2(r,R,t) 



V^7v(r,R,t)/ 



f H{ei) ••• 

H{e2) ■■■ 

■•. : 

; ; ■•. 

V ■■•0 HieN)J 



\ /^i(r,R,t) \ 
^2(r,R,t) 



\^N{r,R,t)J 



(62) 
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where \I/fc is the normed wave function of the k quantum state and 



Hisk) 



h' 



2M 



-A 



R 



r(i-gfc) 

2/i 



-Ar + Vir) 



(63) 



£fc=(^fe|£(|Fi2(r)|)|^fc) 



(64) 



We recall that the average value of any function g of the operators fi, f2, pi, and p2 has 
the following form: 



N 

J ^^ g^dT = J2 \(^kf J % 9[h,r2,M^k),P2{ek)] "^kdr, 
fc=i 



(65) 



where the matrix g is diagonal 



fgiei) 

9(62) 

■•• 

i ; 

V ■■■ 



\ 








(66) 



g{eN)J 



and g{ek) = g[ri,r2,pi{ek),P2{£k)]- In (|65D, the summation is carried on over all discrete 
states of the system. If a continuous spectrum is available, we perform the summation over 
the entire discrete spectrum and the corresponding integration over the whole continuous 
spectrum in (|65|). 

If all Ek tend to zero, all H{ek) become identical and equal to the Hamiltonian of the 
Schrodinger equation, and all wave functions are the eigenfunctions of this Hamiltonian. 



III. DISCRETE SPECTRUM OF A HYDROGENLIKE ATOM WITHIN THE 

MODEL OF NONCOMMUTING OPERATORS OF COORDINATES AND 

MOMENTA OF DIFFERENT PARTICLES 



Consider the discrete spectrum of a hydrogenlike atom. Let two particles with masses mi 
(electron) and m2 (atomic nucleus) be bound by the Coulomb potential V{r) = —Ze'^r~^, 
where Z is the charge of the atomic nucleus. The self-consistent system of integro-differential 
equations for every state with the binding energy Eni (p8D-(p^) can be written as 



h'n 



Snl 



2/i 



dr"^ 



1(1 + iy 



Ze' 



Xniir) = EniXniir), 



(67) 



£nl 



Ze^ 
Ze^ {ii(f 



-dr. 



(68) 



+ 
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xli{r)dr = 1. 



(69) 



Two equations (p7\} and ( |69D are the equations for the normed radial functions of a 
hydrogenhke atom by the Schrodinger theory. Their solutions for bound states are well- 
known (see, e.g., 0): 



Xni{r) = N^y+'F -n + / + 1, 2/ + 2, 



2Zr 



-Zr 



Enirnao^ 



exp 



SniYnaQ^ 



where 



K 



nl 



(2/ + 1)! 



{n + l)\ 



2n{n-l-l)\ 



1/2 



2Z 



«+3/2 



'1 -EniYnaQ^ 



(70) 



(71) 



Here, a,, = h'^f/je^ is the Bohr radius, and F is a degenerate hypergeonietric function. 
Eigenenergies are 



F, 



nl 



fic'^ {aZf 



n^ 



il-£niy 



(72) 



Here, a = e^/hc is the fine structure constant, / = 0, 1, . . . , n — 1, and n = 1, 2, . . . , oo. By 
substituting Xni{r) into Eq. (p5|), we obtain the nonhnear equation for the determination of 

£nl '■ 



Vnl 



Snl I X 

/o 



21+2 



exp (-x) F^ (-n + / + 1, 2/ + 2, x) 1 + 



AiaZf^ 






(73) 



where Vm = I - Sni and Sni = [{21 + 1)!]"' [2n{n - / - 1)!]"' [{n + /)!] 
For a hydrogenlike atom in the ground state, we have: 



Tlu) = - X exp(^-xj 
2 Jo 



4:{aZf 



dx . 



(74) 



With respect to r/iQ, this equation has solutions if aZ < aZ^ = 0.510107 (Fig. |T|). From two 
solutions, a solution being closer to unit is suitable. The second should be omitted since it 
corresponds to the case where the binding energy tends to minus infinity as the interaction 
constant aZ tends to zero, which is physically unacceptable. For Z > Zc, Eq. ([7^ ) has no 
solutions, which means that the given bound state does not exist. 

Fig. displays the binding energy of the ground state of a hydrogenlike atom along with 
the analogous dependence by the Schrodinger theory. For the sake of comparison, we also 
present the corresponding solution of the relativistic Klein-Gordon equation for a spinless 
particle with mass fi in the Coulomb field V{r) = —Ze^r"^ as 



Enl = fiC^ { -1 + 



l + a^Z^ 



n 



2ry2 



/-0.5+ (/ + 0.5)^-a^Z 



1/2 



^1/2^ 



(75) 



We see that the energy levels are situated below the Schrodinger ones and higher than 
that computed according to the Klein-Gordon theory. Analogous calculations can be easily 
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performed for excited states of a hydrogenlike atom (Figs. ^^. In this case, it turns out 
that the Schrodinger levels with given n are split into n closely positioned sublevels since 
the orbital number / can take n values (/ = 0, 1, . . . , n — 1), i.e., the degeneracy is removed. 
All levels with a given n and different I are situated below the corresponding Schrodinger 
level. 

For a hydrogen atom, the parameter of noncommutativity Sni is significantly less than 
1 (£^0 = 0.776 X 10-^£2o = 0.970 x 10-^e2l = 0.324 x 10"^). The splitting of the level 
with n = 2 is also very small and is by two orders of magnitude less than that by the Dirac 
theory for a hydrogen atom. 

The parameter of noncommutativity e of the operators of coordinates and momenta of 
different particles , presented in Fig. ^ decreases as the quantum numbers n and / increase 
(for the same Z), i.e., fully nonrelativistic solutions pass into solutions of the Schrodinger 
equation for large quantum numbers. 

A characteristic feature of the fully nonrelativistic equation is the presence of the critical 
value of the parameter aZ,. for each energy level, which is lacking in the case of the nonrel- 
ativistic Schrodinger equation. For example, for levels with n = 2, aZ^ = 1.401098 for / = 
and aZc = 1.221611 for / = 1. 

As the parameter aZ grows, the average distance between particles decreases. For the 
ground state of a hydrogenlike atom, (|f2 — fi|) = 3^(1 — eioY {2 ficaZ)~^ ^^d takes the 
minimum value (|f2 ~ fi|) ^ 1.33h/fic for Z = Z^. With a further increase in the parameter 



aZ, the self-consistent system of equations (|67D-(|69D has no solutions in the state with n = 1 



i.e., the IS-state cannot exist, and the ground state is the state with n = 2 and I = (2S- 
state) if 0.510 < aZ < 0.847 or 1.222 < aZ < 1.401. In the region 0.847 < aZ < 1.222, the 
ground state is the state with n = 2 and / = 1 since £21 > ^20 (2P-state). Here, we observe the 
substantial distinction from solutions of the nonrelativistic Schrodinger equation, for which 
the ground state is, as known, the IS-state. Solutions of the proposed fully nonrelativistic 
solution for a hydrogenlike atom are also significantly different from solutions of the Klein- 
Gordon equation. The latter has the critical value of the interaction constant aZc = 0.5. 
Above it, the continuous energy spectrum appears, and there "occurs" the so-called fall to 
the center. 

Now we write down the values of quantum Poisson brackets proposed by Dirac P): 

{xi, Plx} = I - ^£nl , (76) 

{X2, P2x} = 1- -J^Snl , (77) 

{Xl, P2x} = ^Snl , (78) 

{X2, Plx] = -J^Snl , (79) 

{xi, X2} = 0, (80) 
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{pix,P2x}=0. (81) 

As Eni -^ 0, these brackets are transformed into the classical Poisson brackets. That is, in 
this case, we have a complete analogy between classical mechanics and quantum mechanics. 
Fig. ^ demonstrates that e significantly differs from zero for systems, whose size is about the 
Compton wavelengths of particles composing the system. In this case, there is no similar 
analogy with classical mechanics. 



IV. CONCLUSION 

The Schrodinger equation for a system of interacting particles is not a strictly nonrela- 
tivistic equation because it is grounded on the implicit assumption about finiteness of the 
interaction propagation velocity. The last means that if the commutator of operators of a 
coordinate and the corresponding momentum of a free particle is defined as 

[x, pj = th , (82) 

this commutator for a system of coupled particles has the same value ih . However, in a 
nonrelativistic quantum system during measurement of the coordinate of a particle, a whole 
transferred momentum is distributed over all particles but is not transferred to only the 
measured one. Therefore, in a system of interacting particles, this commutator should have 
the form 

[x, px] = iM , (83) 

where < 5 < 1 . 

The rejection of the implicit assumption on finiteness of the propagation velocity of 
interactions implies the noncommutativity of the operators of coordinates and momenta of 
different particles. But the operators of coordinates of all particles and operators of momenta 
of all particles mutually commute that allows one to use these collections as independent 
variables. 

The derived self-consistent system of integro-differential equations allows one to separate 
the motion of the center of mass of the system, which moves as a free particle. 

The properties of solutions of the proposed system of equations significantly differ from 
those of Schrodinger solutions for systems, whose size is comparable with the Compton 
wavelength of particles. That is, the consideration of noncommutativity of the operators of 
coordinates and momenta of different particles is important for the quantum mechanics of 
atoms with a large charge of nuclei as well as for the phenomena of nuclear physics, for which 
the size of a system is about the Compton wavelength of particles composing the system. 

In conclusion, the author expresses his gratitude to Dr. V. V. Kukhtin and Prof. I. V. 
Simenog for a very useful discussion of certain problems touched in this work. 
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FIGURES 
FIG. 1. Dependence of the right-hand side of (74) on rj for various values of the parameter aZ. 

FIG. 2. Binding energy Eiq of the ground state of a hydrogenUke atom (72) vs the parameter 
aZ. The upper dotted hue corresponds to the Schrodinger theory, Es = —fic'^{aZ)'^/2, and the 
lower dotted line to the corresponding solution of the Klein-Gordon equation (75). 

FIG. 3. Binding energy E20 of a hydrogenlike atom (72) vs the parameter aZ. The upper 
dotted line corresponds to the Schrodinger theory, Es = —fJ,(P{aZ)'^/8, and the lower dotted line 
to the corresponding solution of the Klein-Gordon equation (75). 

FIG. 4. Binding energy £^21 of a hydrogenlike atom (72) vs the parameter aZ. The upper 
dotted line corresponds to the Schrodinger theory, Es = —fJ,(P{aZ)'^/8, and the lower dotted line 
to the corresponding solution of the Klein-Gordon equation (75). 

FIG. 5. Dependence of the parameter of noncommutativity of operators e for the lowest states 
of a hydrogenlike atom vs the interaction constant aZ. 



16 



1 - {xL<(xL 

2 - aZ = aZ 

3 - aZ > aZ 




0-8 ^,0 ^ 



'10 



|au 


^^^^^^^^^^^^ 










0.1 


> 


\X -v 






\ 


0.2 


- 


\\ 






K-G 'i| 


0.3 


1 1 1 1 1 1 1 1 


1 1 



0.0 0.1 0.2 0.3 0.4 0.5 aZ 



"20 



^ 2 


^^^^^^^^ 




■0.1 


- ^K-a 


^^^^> 

^^^>.» 
^^^>^ 


■0.2 


: 


\; x^ 


■0.3 


- 


\ ^'^ 


■0.4 


- 


\ 


■0.5 


1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 aZ 



'21 



|LlC 



■0.2- 



■0.4- 



■0.6- 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 aZ 




0.15 - 



0.10- 



0.05 - 



0.00 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 aZ 



